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We develop a time-dependent nonequilibrium Green function (NEGF) approach to the problem 
of spin pumping by precessing magnetization in one of the ferromagnetic layers within F|I|F mag- 
netic tunnel junctions (MTJs) or F|I|N semi-MTJs in the presence of intrinsic Rashba spin-orbit 
coupling (SOC) at the F|I interface or the extrinsic SOC in the bulk of F layers of finite thickness 
(F-ferromagnet; N-normal metal; Finsulating barrier). To express the time-averaged pumped charge 
current, or the corresponding dc voltage signal in open circuits that was measured in recent exper- 
iments on semi-MTJs [T. Moriyama et al, Phys. Rev. Lett. 100, 067602 (2008)], we construct 
a novel solution for the double-time-Fourier-transformed NEGFs. The two energy arguments of 
NEGFs in this representation are connected by the Floquet theorem describing multiphoton emis- 
sion and absorption processes. Within this fully quantum-mechanical treatment of the conduction 
electrons, we find that: (i) only in the presence of the interfacial Rashba SOC the non-zero dc pump- 
ing voltage in F|I|N semi-MTJ can emerge at the adiabatic level (i.e., proportional to the microwave 
frequency); (ii) a unique signature of this charge pumping phenomenon, where the Rashba SOC 
resides within the precessing F layer thereby participating in the pumping process, is dc pumping 
voltage that changes sign as the function of the precession cone angle; (Hi) unlike conventional spin 
pumping in the absence of SOCs, where one emitted or absorbed microwave photon is sufficient to 
match the exact solution in the frame rotating with the magnetization, the presence of the Rashba 
SOC requires to take into account up to ten photons in order to reach the asymptotic value of 
pumped charge current; (iv) disorder within F|I|F MTJs can enhance the dc pumping voltage in 
the quasiballistic transport regime; (v) the extrinsic SOC in F|I|F MTJs causes spin relaxation and 
eventually the decay of dc pumping voltage which becomes negligible when the ratio of F layer 
thickness to the spin-diffusion length is around five. 



PACS numbers: 85.75.-d, 72.25.Mk, 72.10.Bg, 73.40.-c 

I. INTRODUCTION 

The spin pumping by precessing magnetization is a 
phenomenon where the moving magnetization of a sin- 
gle ferromagnetic layer, driven by microwave radiation 
under the ferromagnetic resonance conditions (FMR), 
emits spin current into adjacent normal metal layers.^ 
The emitted spin current is pure? in the sense that it is 
not accompanied by any net charge flux. This effect is 
termed "pumping" because it happens in the absence of 
any dc bias voltage, and together with closely related adi- 
abatic quantum pumping of charge^ or spin^ observed in 
quantum dots, falls in the category of problems where an 
open quantum system (i.e., a finite many-particle system 
in contact with particle reservoirs) is exposed to time- 
dependent periodic externals fields. 

Since angular momentum loss carried by emitted pure 
spin current adds extrinsic contribution to Gilbert damp- 
ing, spin pumping has initially been observed^ ^ ^ as an 
increased broadening of FMR spectra upon switching 
from a single ferromagnet (F) layer to F|N multilayers 
(N-normal metal). Therefore, it is also an essential in- 
gredient to understand'' critical current switching in ex- 
periment^ on spin-transfer-torque-driven magnetization 
dynamics. In fact, spin pumping can also be viewed 



as the Onsager reciprocal phenomenorP of spin-transfer 
torque (STT) in which spin current of large enough 
density injected into a ferromagnetic layer either switches 
its magnetization from one static configuration to an- 
other or generates a dynamical situation with steady- 
state precessing magnetization. 

Recent vigorous experimental efforts have focused on 
the direct detection of pure spin current generated by co- 
herent macrospin precession in both ferromagnetic met- 
a^l ^ir f fT3l insulator^I^GS by converting it into dc volt- 
age signal. For example, the experimental techniques 
employed for this purpose include the inverse spin Hall 
effect or the second static F layer as detector within a 
magnetic tunnel junction (MTJ),^^ as well as the precess- 
ing F layer itself which can generate voltag^^at the F|N 
interface by detecting the backflow spin current (due to 
spin accumulation^ driven by pumping into the N layers 
thinner than the spin-diffusion length). 

The origin of the voltage signal of spin pumping in 
F|I|F MTJs can be easily understoocP^^^ as a two step 
process: (1) the magnetization dynamics of the left F 
layer pumps pure spin current across the tunnel barrier 
(I-insulator); (2) the pumped spin current is then filtered 
by the analyzing right F layer whose magnetization is 
static. This generates charge current or, equivalently, dc 
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pumping voltage in an open circuit P^H^^ This voltage is 
proportional to the frequency oc huo of microwaves due 
to adiabatic nature of pumping (in the adiabatic regime, 
formally uj ^ since energy of microwave photons is 
smaller than other relevant energy scales in ferromagnetic 
solids). 

On the other hand, surprisingly large voltage signal 
oc huj observed^^ in F|I|N semi-MTJs, which do not con- 
tain the second analyzing F layer, has remained un- 
explained in virtually all recent theorie^^^HI^^of spin 
pumping in MTJs. Some of these theorie^^-^^ actu- 
ally predict tiny voltage signal which being non-adiabatic 
oc (hw)'^ is the second-order effect. 

Also, there exists several orders of magnitude dis- 
crepancy between underestimated voltage signal of spin 
pumping in F|I|F MTJs obtained in the scattering the- 
oryp^ experimental data,!!^ and overestimated voltage 
signal obtained in the rotating frame approach^^ or the 
tunneling Hamiltonian formalism combined with semi- 
classical modeling of the interplay of spin diffusion and 
self-consistent screening around interfaces. This can be 
traced to different device setups where scattering ap- 
proach was applied to MTJs assuming zercP^ spin ac- 
cumulation in the F layers modeled as semi-infinite leads 
(justified through assumption that spin- flip rate in F is 
larg er th an the tunnel rate), while unrealistically large 
bullP^^ or interfacial^ spin accumulation appears in 
the other two approaches. 

The MTJs employed in spin pumping experiment J^^^^ 
contain F layers of nanoscale thickness whose short spin- 
diffusion lengtlPI' can be modeled by sufficiently strong 
extrinsic spin-orbit scattering. Most importantly, the 
very recent experiments^^ have unveiled a possibility of 
strong Rashba spin-orbit coupling (SOC) at the F|I in- 
terface due to structural inversion asymmetry of the mul- 
tilayered device within which such interfaces reside (e.g., 
Rashba SOC was detected in N|F|I multilayers, but not 
in the N|F|N ones^^). 

However, SOCs have been traditionally neglected in 
a variety of approaches t o spin and charge pumping 
by magnetization dynamics. I^E^l^SISSES r^j^^ SOCs in the 

bulk or on the surface of ferromagnetic materials play 
crucial role in other phenomena, such as the anomalous 
Hall effect-^ or the tunneling anisotropic magnetoresis- 
tanc^^^''^(TAMR). Moreover, the study of the interplay 
between SOC and STT has been recently initiated^^^^^^^ 
throug h the oretical proposals^^ and experimental real- 
ization^SHMI that exploit SOCs for STT-driven magneti- 
zation reversal of a single F layer with greatly reduced 
critical current required when compared to traditional 
spin valves or MTJs with two non-collinear magnetiza- 
tions.EHl 

In the "standard modeF'II'^ of spin pumping in mag- 
netic multilayers containing man}d F and N layers, the 
magnitude of pumped spin current by F|N interfaces is 
computed quantum-mechanically via the Brouwer scat- 
tering formula^^ which then serves as the boundary con- 
dition for the spin-diffusion equatiorPl or enters into 



the so-called circuit theor}^^ where device is split into 
nodes of characteristic size smaller than the spin-diffusion 
length. Thus, in these frameworks SOCs enter only phe- 
nomenologically through finite spin-diffusion length (a 
spin can be flipped by SOC and magnetic impurities in 
N or F layers, as well as by magnon scattering in the F 
layers). 

However, these approaches are not applicable to MTJs 
since the spin accumulation is not well-defined in the 
insulating barrier.^ Moreover, even the magnitude of 
pumped current cannot be obtained from the "standard 
model" formula,^^ governed by the interfacial spin-mixing 
conductance^ when strong SOC is present imme- 
diately at the interface which renders g^-^ an ill-defined 
quantity. 

The recent alternative description^^ of spin pumping 
in F|N multilayers, based on nonequilibrium Green func- 
tion (NEGF) expressions for the local spin and charge 
current densities, has encompassed both the earlier con- 
sidered^ nonlocal diffusion of the spin accumulation at 
the F|N interface generated by magnetization precession 
and the effective field described by the "standard model" 
(where spin accumulation does not build at the interface 
since spin-flip relaxation rate is assumed to be sufficiently 
larger than the spin injection rate). However, this frame- 
worl;^ has treated SOCs only in the N layer away from 
the precessing F layer in order to analyze how each of 
these two pumped spin currents can be converted into 
charge current by the inverse spin Hall effect due to the 
extrinsic or intrinsic SOCs. That is, in this theory SOCs 
are not essential for the discussion of spin pumping effect 
itself. 

We note that other groups have also recently identi- 
fied the importance of adding SOCs explicitly into the 
description of spin and charge pumping by magnetiza- 
tion dynamics. For example, a generalized scattering 
theory of adiabatic charge pumping by a single precess- 
ing F layer within N|F|N junctions containing SOCs has 
been formulated in Ref. 9, Also, the pumping of current 
of magnetic monopoles and the associated (via Ampere 
law) charge current flowing in the plane of the Rashba 
SO-coupled interface (rather than perpendicular to it as 
is the focus of our study) between the precessing F layer 
and a nonmagnetic layer has been predicted in Ref. |39l 

Here we develop a NEGF-based quantum transport 
theory of spin current pumping, its propagation, and 
conversion into electrically measurable signals in F|I|F 
and F|I|N junctions depicted in Fig. [l] The junctions 
are described by the microscopic time-dependent (due to 
precessing magnetization of one of the F layers) Hamil- 
tonian which is tailored to take into account nanoscale 
thickness of F layers within which we include terms de- 
scribing disorder and extrinsic SOC in the F layers, as 
well as possibly strong Rashba SOC at the F|I inter- 
faces. Our theory starts from the equations of motion 
generated by such Hamiltonian for NEGFs which de- 
pend on two time variables, and then finds a computa- 
tionally efficient solution to such equations which physi- 
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FIG. 1: (Color online) (a) F|I|F MTJ and (b) F|I|N semi-MTJ 
with precessing magnetization of a single F layer are modeled 
on a simple cubic finite-size tight-binding lattice attached to 
semi-infinite ideal (disorder and interaction free) N leads. The 
thicknesses of the ferromagnetic layers and thin insulating 
barrier is measured using the number of atomic monolayers 
dp and dj, respectively. For example, dp = S and = 4 
in the illustration, while in the actual calculations we use 
c/f = 50 and = 5 monolayers of cross section 20 x 20 lattice 
sites. The interfacial Rashba SOC due to structural inversion 
asymmetry of the junction is included in the last monolayer 
of the F slab that is in direct contact with the tunnel barrier 
I. The F layers can also include disorder modeled as a random 
on-site potential and the corresponding extrinsic SOC, while 
binary alloy disorder in the I layer models AlO^-type tunnel 
barrier. 



cally describes processes where a specific number of mi- 
crowave photons is absorbed or emitted by propagat- 
ing electrons in the course of pumping. This solution 
for time-dependent NEGFs allows us to obtain time- 
averaged total charge current in the N leads of the junc- 
tions shown in Fig. [l] or the corresponding dc pumping 
voltage in the corresponding open circuits. The formulas 
which we derive for pumped currents are also applicable 
to any problem where a quantum-mechanical system is 
exposed to periodic time-dependent external field, inde- 
pendently of its frequency (i.e., including both adiabatic 
and non-adiabatic regimes) or amplitude (i.e., includ- 
ing both perturbative and non-perturbative regimes), as 
long as electron-electron, electron-phonon and electron- 
magnon interactions can be neglected. 

The paper is organized as follows. In Sec. [n| we dis- 
cuss how to tune parameters of the microscopic Hamil- 



tonian in order to reproduce properties of MTJs and 
semi-MTJs employed in experiments,'^^'^ such as their 
tunneling magnetoresistance, diffusive nature of trans- 
port within the F layers, finite spin-diffusion length in 
F layers and strong interfacial Rashba SOC. Section III 
discusses NEGF equations of motion and how to solve 
such equations after converting them into algebraic ones 
via double Fourier transform. In this section we also 
employ clean F|I|F junctions (with semi-infinite or finite 
thickness F layers) as a test bed to compare our theory 
to the scattering formulas for pumping in MTJs derived 
in Ref. |16l as well as to previously developed^^ solution 
to double-time-Fourier-transformed NEGF equations us- 



ing continued fractions. In Sec.[IV|we discuss properties 
of time-averaged pumped charge current and the corre- 
sponding dc voltage signal in open circuits for clean F|I|N 
and F|I|F junctions as a function of the strength of in- 
terfacial Rashba SOC. The effect of disorder in F and 
I layers, as well as the extrinsic SOC within F layers, 
on the dc pumping voltage in F|I|F MTJs is discussed 



in Sec. [Vl We conclude in Sec. [Vl) Readers seeking to 
understand physical effects of SOC on microwave-driven 
MTJs may wish to start with Sec. Inland then jump to 



Sec. IV and Sec. rvl while those interested in theoretical 



II. MTJ DEVICE SETUP AND ITS 
HAMILTONIAN 



The MTJ and semi-MTJ we study are illustrated in 
Fig. [TJa) and (b), respectively. Each atomic monolayer 
shown in Fig. [l] is modeled on the square tight-binding 
lattice with single s-orbital per site. Since disorder is in- 
cluded as random potential in real space, atomic mono- 
layers are of finite size 20 x 20 sites. This can also be 
viewed as the cross section of a supercell which is re- 
peated periodically in the transverse direction,!^ while 
its size is sufficient to allow one to compute all quan- 
tities at the F point (i.e., without the need to perform 
k-point sampling). We have checked that dc pumping 
voltage remains constant as one increases cross section 
size beyond 20 x 20 sites. 

The ferromagnetic layers consist of dp = 50 such 
monolayers, so that their thickness is o:^ 15 nm (assum- 
ing typical lattice spacing a 3 A), which closely mim- 
ics F layers employed experimentally.'^^'^ The thickness 
of the insulating barrier is dj = 5 atomic monolayers. 
The finite-size F|I|F or F|I|N multilayer is connected to 
macroscopic reservoirs via two semi-infinite ideal (i.e., 
disorder, spin and charge interaction- free) N leads to 
form a two-terminal device required for both NEGF and 
scattering theory analysis. 

The general time-dependent Hamiltonian describing 
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these two devices can be written as 

{ii'), ss' ii' ,ss' 

+ E 4stw^i's'- (1) 

{u'),ss' 

Its time dependence stems from the unit vector mi(t) 
along the local magnetization direction within the left 
F layer, which is assumed to be spatially uniform and 
steadily precessing around the z axis with a constant cone 
angle. The value of angle is controlled by the input 
microwave power (typically < 20° in the recent experi- 
mental^. The operators Cj^ (cicr) create (annihilate) elec- 
tron with spin a at site i = {ix^iy^iz)^ and 7 is the nearest 
neighbor hopping which sets the unit energy scale. The 
coupling of itinerant electrons to collective magnetization 
dynamics is described through the material-dependent 
exchange potential Ai, where & = {cFx^ o-y^&z) is the vec- 
tor of the Pauli matrices and [d-]ss' denotes the Pauli 
matrix elements. 

The disorder within F layers can be introduced 
using the uniformly distributed random variable 
s[ G [—W/2^W/2] which models isotropic short-range 
static impurity potential. To account for the proper- 
ties of amorphous AlO^c tunnel barrier in MTJs and 
semi-MTJs employed in the recent spin pumping experi- 
mentspl^the on-site potential on I monolayers is chosen 
as = Ui) ± 6Ui) where random fluctuations SUi) mimic 
binary alloy disorder. 

The impurity potential in the F layers also generates 
extrinsic SOC, as described by the third sum in Eq. 0. 
This can be viewed as the lattice version^^ of the Thomas 
term A(t7 x p) • VVdis in the Pauli-Schrodinger equation 
so that 

That is, on the tight-binding lattice the extrinsic SOC 
acts as additional spin-dependent hopping between both 
nearest neighbor and next nearest neighbor sites. Here 
{ex^ey^e^) are the unit vectors along the x^y^z axis re- 
spectively and Aeso = A/4a. 

The fourth sum in Eq. ([T]) is the tight-binding rep- 
resentation^^ of the Rashba SOC written in terms of a 
generalized nearest neighbor hopping term that acts as 
2x2 Hermitian matrix in the spin space: 

4. _ / -^Rso^z (i = i' + e^) . . 

"'"l+mso^, (i = i' + e,) ' 

The continuous version of the Rashba SOC, 
<^RSo(<^ X p) • e^^/fi, has been traditionally studied 
in the context of two-dimensional electron gases 
(2DEGs) within semiconductor heterostructures with 




Disorder Strength W/y 

FIG. 2: (Color online) The Fano factor of the zero- 
temperature and zero-frequency shot noise vs. the disorder 
strength W in transport through N slab with static disorder 
or F slab with both static disorder and the corresponding ex- 
trinsic SOC of strength Aeso = O.O47. Both slabs consists of 
50 monolayers (containing 20 x 20 atoms per cross section) 
which are connected to two semi-infinite ideal N leads. 



structural inversion asymmetry in the growth direc- 
tion.^^ Nevertheless, several experiments have recently 
reported evidence of the Rashba SOC-induced split- 
ting of the surface states in both non-magnetic and 
magnetic metals using angle-resolved photoelectron 
spectroscopy.^^ The very recent transport experiments^ 
have demonstrated Rashba SOC-induced STT of a 
single thin F layer embedded between two asymmetric 
interfaces. For example, such effect was observed in 
Pt|Co|A10a; multilayers, but not in the inversion sym- 
metric ones Ft I Co I Ft. The experiment of Ref. 22 has 
also utilized heavy atoms and surface oxidation to create 
strong out-of-plane potential gradient in Pt|Co|A10aj 
junctions and enhance the interfacial Rashba SOC. 

This motivates the introduction of the Rashba SOC 
term of strength 7rso into the Hamiltonian in Eq. ([T]), 
which we set to be non-zero only on the last monolayer of 
the precessing F layer that is in the direct contact with 
the first monolayer of the tunnel barrier, as illustrated 
in Fig. [1] Note that the exact location of the interfacial 
Rashba SOC eventually requires fitting the Hamiltonian 
parameters to first-principles analysis. ^^-^ 



A. How to tune Hamiltonian parameters to 
reproduce properties of MTJs used in experiments 

The diffusive transport regime within F or N layers 
is defined semiclassically by the requirement that mean 
free path £ is smaller then the thickness of such layers, 
£ <^ dp OT £ <^ In quantum transport calculations, 
the easiest way to select proper range of values for the 
disorder strength W which ensures diffusive regime is to 
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compute the Fano factor F = S/2eI of the shot noise 
whose zero-temperature and zero-frequency noise power 
is S. For the diffusive metahic wires, F = 1/3 is universal 
in the sense of being independent of the impurity distri- 
bution, band structure, and shape of the conductor. In 
Fig. |2] we plot the Fano factor with increasing disorder 
strength W for both conventional disordered N layer at- 
tached to two ideal semi-infinite leads and F layer of the 
same dimensions with both disorder and extrinsic SOC. 
Using Fig. [2j we select W = to ensure semiclassical 
diffusive transport regime. 

Note that spin-dependent effects on the shot noise 
are reveled only when spin-polarized current is injected 
and the corresponding spin-resolved Fano factors are de- 
fined.'^ Thus, for unpolarized injected charge current, 
both F and N layers have virtually the same Fano factor 
in the diffusive regime, while larger Fano factor for the 
F layer in the quasiballistic regime is due to increased 
scattering at the N|F interfaces because of non-zero A. 



1. Extrinsic SOC in the bulk of F layers 

In both N and F layers, spin-fiip scattering will destroy 
nonequilibrium spin accumulation which is typically 
accounted^ through phenomenological spin-diffusion 
length I/gf. Over this length scale, an injected spin ac- 
cumulation loses its polarization so that Lgf in ferromag- 
nets defines the magnetoelectrically active region of F 
layer in contact with N layer. In metallic ferromagnetic 
materials, Lgf rangeJi^ from 5 nm in Ni8oFe2o (permal- 
loy) , which is often employed as precessing F layer in spin 
pumping experiments,^^ to 50 nm in Co. 

To understand how to tune the strength Aeso of the 
extrinsic SOC term in Eq. ([T]) in order to generate dif- 
ferent experimental situations in F layers of nanoscale 
thickness, such as Lgf < dp or Lgf > dp (where spin-fiip 
processes essentially become unimportant), we compute 
the spin density matrix^^ of collected transported spins in 
the right lead 2 after fully spin-polarized charge current 
is injected from the left lead 1 across the F layer: 




The elements of the transmission matrix t^f in this 
formula determine the probability |[t2f ]nmP for spin- 
cr' electron incident in lead 1 in the orbital conducting 
channel |m) to be transmitted to lead 2 as spin-a elec- 
tron in channel |n). Therefore, such amplitudes also de- 
termine the corresponding spin-resolved conductances, 

G21 = T^n,m=i IMiUml^- The spiu density matrix 
makes it possible to extract the detected current polariza- 
tion vector P°^* = (P°^SP°^SP°^*) whose magnitude 
pout _ I pout I g^ygg ^i^Q so-called current polarization 
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FIG. 3: The decay of current polarization along the diffusive 
F layer with static disorder of strength W = 3j and the ex- 
trinsic SOC of strength Aeso = O.I7. The F layer is attached 
to two semi- infinite ideal N leads where charge current which 
is 100 % spin-polarized along the z-axis, P^^ = (0, 0, 1), is in- 
jected from the left N lead and P°^* is computed in the right 
N lead for F layer of thickness dp- The unit vector of the 
magnetization in F is either parallel (solid line) or antiparal- 
lel (dashed line) to the z-axis. The inset shows spin-diffusion 
length as a function of Aeso when F layer is replaced by diffu- 
sive N layer with different strengths of extrinsic SOC, where 
each value of Lsf is extracted by fitting exponentially decaying 
function to P°^^ vs. dN curves. 



measured experimentally.^^ The matrices and 
are computed for the two-terminal N|F|N device where F 
layer with static magnetization pointing along the z-axis 
and disorder of strength = 87 is embedded between 
two semi-infinite N leads [this means that the number of 
conducting channels M in Eq. Q is M = 400 for the 
selected cross section of 20 x 20 sites]. 

The measured current polarization of permalloy at 
room temperature ranges from P = 0.32 to P = 0.5, de- 
pending on the experimental technique employed. Since 
we find only tiny fluctuations of P^;, Py ^ 10~^ in the 
presence of non-zero A and Aeso, we use P°^* = |P°^*| 
as the measure of current polarization. We first tune 
A = 27 of the F layer with Aeso = to obtain P°^* 0.5 
at the Fermi energy Ef = —37. Then we compute the de- 
cay of P°^* with increasing length of the diffusive F layer 
with non-zero Aeso, as shown in Fig. |3] For sufficiently 
thick F layer, these curves saturate at |P°^*| — 0.5. On 
the other hand, the same calculation for the diffusive N 
layer with non-zero Aeso gives usual exponentially de- 
caying P°^* vs. dN curves due to spin diffusion, whose 
fitting establishes the correspondence between Aeso val- 
ues used in our study and microscopically determined 
spin-diffusion length Lgf. The dependence Pgf oc 1/Aeso 
shown in the inset in Fig. |3] is expected for the diffusive 
transport regime. 

With A and Ef specified in this fashion to ensure that 
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FIG. 4: (Color online) (a) The setup for the measurement 
of the out-of-plane TAMR in F|I|N semi-MTJ, defined by 
Eq. (|6| , as a function of the angle between the static magne- 
tization of the F layer and the transport direction (the x-axis). 
In panel (b), the Rashba SOC at the F monolayer in contact 
with the tunnel barrier I is fixed at 7rso = O.57, while panel 
(c) shows TAMR((/) = 90°) for different values of 7rso. (d) 
The setup for the measurement of the out-of-plane TAMR in 
F|I|F MTJ, defined by Eq. ([7|, as a function of the magne- 
tization orientation in each of the two F layers with respect 
to the transport direction. The TAMR depends on the abso- 
lute magnetization directions mi and m2. In panel (e), the 
Rashba SOC of strength 7rso = 0.67 is present at both F 
monolayers in contact with the tunnel barrier I. The F layers 
in both semi-MTJ and MTJ have finite thickness dp — 50. 



current polarization of permalloy slab matches experi- 
mentally measured values, we finally tune the height of 
the potential barrier [/^ = 97 in the tunnel barrier I of 
thickness dj = 5 to tune the "optimistic" tunneling mag- 
netoresistance (TMR) defined as 



TMR 



Rap — Rp 
Rp ' 



(5) 



to TMR=50 % for the F|I|F MTJs employed experimen- 
tallyP^^ Here Rap is the resistance with antiparallel 
configuration of magnetizations in the F layers of thick- 
ness dp = 50, while Rp is the junction resistance when 
magnetizations are parallel. Since both of these resis- 
tances are dominated by the tunnel barrier potential, 
they are computed for clean j unctions. To model AlO^c 
tunnel barrier, we use binary alloy disorder character- 
izecPlby SUb = 0.5j. 



2. TAMR and spin dephasing in perpendicular transport 
through interfaces with the Rashba SOC 

To understand the correspondence between the 
strength of the Rashba SOC measured by the spin- 
dependent hopping parameter 7rso in Eq. (|3| and the 
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FIG. 5: (Color online) The magnitude |P°^'| of the spin- 
polarization vector of outgoing charge current in the right N 
lead after fully |P^^| = 1 spin-polarized current is injected 
from the left N lead traversing a monolayer with the Rashba 
SOC of strength 7rso. The spin-polarization vector P^^ can 
point along three different axes of the coordinate system in 
Fig. [1] where P^"" = (0,1,0) and P^"" = (0,0,1) are parallel 
to the Rashba monolayer while P^^ = (1,0,0) is orthogo- 
nal to the Rashba monolayer. The direction of p°^* remains 
collinear with P^^, as illustrated in the inset. 



values encountered in experimental devices,'^ we com- 
pute the so-called out-of-plane TAMR coefficient for 
FIIIN semi-MTJ which is defined a^ 



TAMR (0) 



R{(l)) - R{0) 



(6) 



for the device setup illustrated in Fig. [4]^a). In dc trans- 
port measurements of TAMR, the magnetization direc- 
tion in F layer provides a control knob orienting the spin, 
while the magnetic anisotropy is determined by the inter- 
face symmetry rather than by the symmetry of the bulk 
materials. Here i^(0) is the resistance of semi-MTJ when 
static magnetization of its F layer is parallel to the x-axis 
as the direction of transport in Fig.|4ja), and R{(t)) is the 
junction resistance when magnetization is rotated by an 
angle (j) with respect to the x-axis within the xz-plane. 

Figure gb) shows TAMR(0) at fixed Rashba SOC, 
while the maximum TAMR(0 = 90°) vs. the strength 
of Rashba SOC is plotted in Fig. ^c). Compared to the 
weak Rashba SOC in 2DEGs where typically 7rso — 
O.OI7, the interfacial Rashba SOC in semi-MTJs has to 
be rater strong (as achieved in the recent experiments^ 
in order to generate observable TAMR. Since the inter- 
facial SOI is linear in momentum, TAMR vanishes at the 
first order in 7rso after averaging over the Fermi sphere. 
However the ferromagnet contains local exchange field 
and a net transfer of angular momentum occurs at the 
second order, so that TAMR oc 7^30- This region of 
small TAMR occurs for 7rso ^ O.47 in Fig.jlj^c), beyond 
which higher order processes start to play the role and 
TAMR increases faster with increasing 7rso- 
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For F|I|F junctions, one can define the out-of-plane 
TAMR coefficient a&^i 



current by a remarkably simple formula 



TAMR (6*, (/)) = 



R{e,<p)- R{e,{)) 
R{e,Q) ' 



(7) 



where the meaning of angles 6 and (j) is explained in 
Fig. ^d). Since TAMR coefficient for F|I|N has only 
one angle argument, there is no ambiguity in using the 
same TAMR notation for both cases. The out-of-plane 
TAMR for F|I|F MTJ is shown in Fig. ge). 

Unli ke a mply studied lateral spin transport in 
2DEGJ^^'^ or interfaced under the influence of the 
Rashba SOCj^ very little is known about the effect of 
such interfaces on spin transport perpendicular to the 
plane, as illustrated by the measurement geometry in the 
inset of Fig. [sj We clarify their effect by using Eq. ^ to 
obtain the spin-polarization vector P^^* of the current in 
the right N lead after 100% spin-polarized charge current 
with |P^"| = 1 is injected from the left N lead. The result 
in Fig. [5] shows spin dephasing, where the outgoing spin 
polarization vector po^* remains in the same direction as 
P^", but with reduced magnitude |po^*| < 1. The degree 
of dephasing depends on the direction (perpendicular or 
parallel) of the initial spin polarization with respect to 
the Rashba interface. 



III. NEGF APPROACH TO PUMPING BY 
PRECESSING MAGNETIZATION 

Theoretical studies of quantum charge pumping in 
noninteracting phase-coherent systems have been con- 
ducted using a variety of approaches. In the adiabatic 
regime, Brouwer scattering formulsP^ is often used as an 
elegant geometrical description of the charge pumped per 
cycle in terms of the instantaneous scattering matrices 
of the system. The adiabatic regime occurs when time- 
dependence of the driving field parameters is slow in com- 
parison to the characteristic time scales of the system, 
such as the electron dwell time, so that electrons traverse 
the device as if the external potential landscape if frozen 
in time. Approaches beyond adiabatic regime include 
Floquet scattering theory,!^ iterative solutions of time- 
dependent states^ and variations of the NEGF formal- 
ism.SSI^ Moreover, the generality of the time-dependent 
NEGF framework^^ makes it a usual choice in the studies 
of pumping in the presence of strong Coulomb interac- 
tions.!^ 

Among these approaches, Brouwer scattering formular^ 
and NEGF formalism^^ have been employed to describe 
experiments on spin pumping by moving magnetization 
in magnetic multilayers. Unlike quantum charge pump- 
ing, the spin pumping in magnetic multilayers is robust 
and ubiquitous effect at room temperature. Nevertheless, 
the match between Brouwer scattering formulsfll and ex- 
periments on F|N multilayers is excellent due to pumped 
spin current being determined by the processes at the 
F|N interface. The scattering theory expresses pumped 



j5p ^ 



47r 



Re^ 



dm 

lit' 



whose dc component is given by 



hjuj 



g ' ^ sm 



(8) 



(9) 



However, the derivation leading to this formula, as well 
as the very definition of the spin-mixing conductance^ 
g^^ (where Re is its real part) of the F|N interface, 
assumes absence of any spin- flips. ^ 

Naively, one could numerically evaluate the Brouwer 
scattering formula for the whole device without intro- 
ducing ill-defined g'^^ in the presence of interfacial SOCs. 
However, SOC renders all components of pumped spin 
or charge current time-dependent so that one has to 
compute the scattering matrix at all times within one 
period r = 27r/a; of the pumping cycle an d th en find 
the time-averaged value of pumped currents.l^^ This is 
prohibitively expensive for 3D system composed of large 
number of atomic orbitals [such as the device in Fig. [TJa) 
whose Hamiltonian matrix is of the size 84000 x 84000], 
especially in the presence of disorder where additional 
averaging over impurity configurations is required. Al- 
though this could be achieved for smaller device sizes, 
we find that the maximum value of pumped current os- 
cillating in time is orders of magnitude larger than its 
average value over one period in the case of MTJs which 
prevents the estimate of experimentally relevant time- 
averaged values from the numerical data. 

The rotating frame approach j ^^*^^ where pumping due 
to precessing magnetization is mapped onto a dc trans- 
port within a four-terminal device whose currents can 
be computecF^ using NEGFs, is also inapplicable in the 
presence of SOCs or other spin-flip mechanisms. This 
is due to the fact that the same unitary transformation 
(discussed in Sec. [HI B ) which maps time-dependent Zee- 
man term in Eq. Q to the one frozen at t = 0, generates 
new time-dependent SOC terms in the rotating frame. 

The time-dependent NEGF formalismP^^ makes it 
possible to obtain directly the time- averaged current. 
The intricacy in solving its equations stems from the fact 
that nonequilibrium problems are not time-translation 
invariant, so one has to work with equations for most 
general case of NEGFs which depend on two time vari- 
ables. There are two independent GFs that need to be 
determined: (1) the retarded GF 

G\i',sAt,t') = -ie{t-t'){{cUt),clAt')}), (10) 

describes the density of available quantum-mechanical 
states; and (2) the lesser GF 



(11) 



determines how electrons occupy those quantum states. 
Here (...) denotes the nonequilibrium statistical aver- 
age, and we use h = 1 to simplify notation in this 
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Section. Both GFs can be extracted from the contour- 
ordered GF defined for any two-time values that lie on 
the Keldysh contoui!^ (consisting of a two-way path that 
begins at — oo time, draws forward to +00, and then 



backward to —00). 

The retarded GF is governed by the following equation 
of motioi]^ 



+00 

{t,t')-U{t)G''{t,t')- J dt"i:'-{t-t")G''{t",t') = S{t-t'), 



(12) 



where we use G^, H(t), and notation to emphasize that these are matrices whose indices represent space and 

spin degrees of freedom. In noninteracting systems, the retarded self-energy T,^{t) = ^p{t) is simply the sum of 
self-energies due to leads p attached to the sample. 

The lesser GF satisfies the Keldysh integral equation 



+ CO +CO 



'~it,t')= J dh j dt2G'^(t,ti)S<(ti,t2)G"(t2,0, 



(13) 



where the advanced GF is related to the retarded one through G^(t,t') = [G^{t' ^t)\^ . In the case of noninteracting 
systems, the lesser self-energy = ^p^pit) is only due to attached leads, which can be obtained from the 

retarded ones using X)^ = ifpTp(t). Here fp is the Fermi function of the macroscopic reservoir to which the lead p is 
attached at infinity and Vp(t) = i{^p — [T,p{t)]^). 

The Heisenberg equation for the charge Q = e^j^^ Cj^Cig/ or spin density Sa = l/2XliGp ss' ^lsW(^]ss'Cis' opera- 
tors of electrons in lead p then yields expressions for time-dependent total charge current^^ 



/ . \ +00 
Ipit) = = -2eRe J dt'Tr [G'-(t,OS<(t' - t) + G<(t,t')S°(t' - t)], 



(14) 



or spin current 




+00 



(15) 



in lead p. Here T,^{t) = [5];(-t)]t. The integration of 
these expressions, limT^oo ^ JIt%^''')^ yields the cor- 
responding time-averaged quantities Ip and . 



In stationary problems G^(t,t') and G^(t,t') depend 
on the time-difference t—t'^ which allows to Fourier trans- 
form them into functions of a single energy argument and 
reduce the set of coupled integral and integro-differential 
equations to a set of algebraic equations. On the other 
hand, when the device Hamiltonian depends on time ex- 
plicitly, one has to work with both times. Since directly 
solving equations Eq. (12) and (13) is cumbersome, it is 



advantageous to switch to a more convenient representa- 
tion. The typical choices used for problems containing 
periodic time-dependent fields are: (i) the double-time 



Fourier transfornPEZEHl 



G^'<(t,tO 



+00 +00 
dE f E' 



27r / 27r 



(16) 

{ii) the single Fourier transfomi^^ in the time difference 
t-t' 



+ CX) 



G'''<(t,t') 



— e-'^(*-*')G'-'<(t,^); (17) 
27r 



and (m) the so-called Floquet matrix form GJ^^(co') de- 
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fined hy^ 



+ 00 ^ 



/2 
/2 

m-\-n 



(18) 



The expressions in Eq. (18), where t^Qi = t — t' and 



^av = (^ + ^0/2? exploit the periodicity condition G(t + 
r, + r) = G{t^t') and are, therefore, GF counter- 
part of the Floquet matrix representation for periodicahy 
time-dependent Hamiltonian H(t + r) = H{t) and its 
eigenstates (with the "Brihouin zone" of energies being 
-Vt/2 <E < Vt/2). 



A. Exact multiphoton solution to 
double-time- Fourier-transformed NEGFs 



Here we adopt the double-time Fourier transform in 
Eq. (16), which has been used frequently to solve NEGF 
equatio ns for non-adiabatic charge pumpin^l or spin 
pumpin#3^ from the 2DEG with the Rashba SOC 
driven by time-periodic external fields. Due to the Flo- 
quet theorem, the double-time-Fourier-transformed re- 
tarded GF G'^{E,E') must take the form 



G\E, E') = G\E, E^nuj) = G;(E). (19) 

The coupling of energies E and E -\- nu (n is integer) 
indicates how multiphoton exchange processes contribute 
toward the pumped current. 



The double-time- Fourier-transformed Eq. ( 12 ) is given 



by: 



+ 00 



EG'^iE.E 



dE' 



Il{E-E')G'^{E\E^nuj) 



(20) 



-5]^(£;)G^(£;, E^nLj)= 27T5{nuj) 

The Fourier transform of the Hamiltonian 

ll{E) = 27r[lloS{E) + \S{E + cj) + V^J(£; - a;)], (21) 

consists of the first term due to the time-independent 
part Ho of Eq. ([T]) , while the other two terms are Fourier- 
transformed harmonic time-dependent part cast in the 
form H'(t) = Ve'^* + V^-'^K Here V is the matrix 
representation of the operator 

V = - ^^Ai sine {[a^jss' - i[cry]ss') cls^is' ^ (22) 

extracted from Eq. ([T]) as the term carrying the periodic 
time-dependence. 



By substituting Eq. ([21| into Eq. ( |2Q| ) we arrive at the 
following equation: 



[Ei^n-ii- t'^iE + n)]G'^{E) = i. 



(23) 



To simplify the notation, we use 



H 




Ho V 
Vt Ho V 

Vt Ho 





(24) 



and 



n = 



2cj1 

-ujl 



ul 

2cj1 



(25) 



Here symbol A is used to denote a matrix which acts in 
the Hilbert space Hei Hph, where the dimension of the 
Hilbert space of photons Hph is infinite. The unit matrix 
in the Hilbert space of a single electron Hei is 1, and the 
unit matrix in Hei Uph is denoted by i. 

Since higher order multiphoton processes yield pro- 
gressively smaller contribution to the pumped current, 
we restrict the dimension of Hph by considering up to 
A/'ph photons. In this case, the dimensionality of Hph is 
2A/ph + 1 since one can have processes with no photon 
exchange n = or the maximum of n = A^ph photons is 
absorbed or emitted. This means that if we keep only 
^ph = 1 processes, G{E) is a matrix of the dimension 
/sites X 2^ X 3ph where /gites is the number of lattice sites 
within the central region of devices in Fig.[l] 2^ takes into 
account the spin degrees of freedom, and 3ph is due to 
single (or no) photon exchange in the course of pumping. 

The Keldysh equation (13) in this representation is 
written as: 



G^{E) = G'^{E)t' 



{E^n)G\E). 



(26) 



The knowledge of G^{E) makes it possible to obtain the 
time-averaged total pumped charge current in lead p in 
the absence of any dc bias voltage 



ph 



+00 

J dETr[tpiG''t G" 



fpG'^ffG"], (27) 



^ Tp and f 



f{E + f2). Since the trace in 
iTrph, is summing over contribu- 



where f = ^ 
the integrand, Tr = Tr, 
tions from different photon exchange processes, the de- 
nominator includes 2A^ph to avoid double counting. Note 
that the part of the trace operating in Hph space en- 
sures the current conservation in our solution to NEGF 
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equations. The analogous formula for the pumped spin 
current into lead p in the absence of any dc bias voltage 
is given by 



1 



+00 



4K 



ph 



dETi [aa{tptG''t G"" - fpG f fG")j;28) 



- r - ^ - as 



Equations (27) and (28) are the central outcome of our 



formalism and can be applied to arbitrary charge or spin 
pumping problem. 

For the specific problem of pumping by precessing mag- 
netization driven by microwaves, we take into account 
that hu <^ Ef and simplify Eq. (27) accordingly by ex- 
panding the Fermi function 



f = f{EF)i-mf/dE. 



(29) 



This leads to the following adiabatic expression for 
pumped charge current: 



2K 



-TT{tp{EF)[n,G\EF)t{EF)]G"{EF)} 



ph 



assuming zero temperature. The commutator 



(30) 



[Q,G {EF)t{EF)] = ^G {EF)r{EF)-G {EF)t{EF)^ 

(31) 

allows us to make the notation more compact. 

The time-averaged value of the pumped charge cur- 
rent is translated into the dc pumping voltage in an open 
circuit via 



^ump — 



G{oy 



(32) 



which is th e quantity measured in the recent experi- 
ments.E^Here G(l9) is the conductance of F|I|F or F|I|N 
junctions computed by tilting the static magnetization of 
the first F layer by an angle away from the z-axis and 
by applying the linear-response bias voltage between the 
N leads attached to the junction. 



B. Comparison with the rotating frame approach 

In the absence of interactions causing spin-flips, such as 
the SOC, it is possible to convert the complicated time- 
dependent transport problem posed by the presence of 
precessing magnetization in the devices in Fig. [ Tjinto the 
time-independent one by performing the unitar y^^^-*^ * ^^ * 
transformation of the Hamiltonian Eq. ([T]) 

i^rot = UH{t)U^ - ^^U^U^ = H(t = 0) - |(7,. (33) 

Here the unitary operator is given hy U = g*a;cr^^/2 
for m(t) precessing counterclockwise. The transformed 
Hamiltonian Hrot is time-independent in the frame ro- 
tating with the magnetization. The Zeeman term 



which emerges uniformly in the central region of devices 
in Fig.[T]and their N leads, will spin-split the bands of the 
N leads, thereby providing an intuitively appealing rotat- 
ing frame^^ picture of pumping. In this picture, the N 
leads in the rotating frame are labeled by (p, a) [p = L^R 
and a =t,i] and they are biased by the electrochemical 
potential differences ii^j^—ji^^ = ou and /i^— /i^ = uj. Thus, 
these leads behave as effective half-metallic ferromagnets 
which emit or absorb only one spin species. The counter- 
propagating dc currents of spin-polarized electrons flow- 
ing from lead /i^ to lead /i^,, where electrons precess in 

the magnetic fleld of H{t = 0) frozen at an angle with 
respect to the 2;-axis in order to enter into oppositely po- 
larized lead, can be computed using NEGF,^^ transmis- 
sion matrices, or the tunneling Hamiltonian approach.C3 
However, the rotating frame approach cannot be ap- 
plied to systems containing SOCs (or any other source 
of spin flips) because unitary transformation would gen- 
erate time-dependent SOC terms in Eq. (33). Never- 



theless, it serves as a useful tool to compare the range 
of validity of different pumping formulas because the 
transport problem deflned by Hamiltonian Eq. ([T]) with 
7rso = ^ESO = can be solved exactly in the rotat- 
ing frame. We start by flrst extending the NEGF-based 
formulas for pumped current in the rotating frame for 
devices with semi-inflnite N leads treated in Ref. [18] to 
those with semi-inflnite F leads. This eventually makes 
it possible to understand the origin of the orders of mag- 
nitude discrepancy between predictions made in Ref. |16] 
and Ref.ilSjfor the dc pumping voltages in the same type 
of MTJs. 

Since the system in the rotating frame is stationary, 
NEGFs which depend there only on the time difference 
t — t' can be Fourier transformed to work with functions 
of a single energy argument 



[^l-H,ot-SLt(^)]" 



{E)GUE)- 



(34) 
(35) 



The retarded self energy in the rotating frame is obtained 
from the self-energy in the lab frame through a simple 
shift of its argument 

i:U^^{E) = ^(^E+^-coa,y (36) 

The lesser self-energy in the rotating frame is then given 

by 



^<,^p{E) = if(^E+ ^cva,^ Fp (^E + \u 



(37) 



This leads to G\^^{E) = G''{E + uja^/2) for the retarded 
GF and G<,(£:) = i G:,,{rp{E + uaj2)G:,, for the 
lesser one, where f = f {E -\-uj 0-^/2) is the Fermi function 
in the rotating frame written as a 2 x 2 matrix in the 
spin space. Using GJ^q^(^), GJ^^(^) and the following 
identity 



^rot 
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leads to an expression for pumped charge current in lead 
P 



+00 



(39) 



where T = ^^Tp. Here = f{E + and = 

f{E — uj/2) are the diagonal elements of i{E). 

Thus, according to Eq. (39) only electrons whose en- 
ergies fall into the interval [Ef — ^7 + ou] participate 
in pumping (at zero temperature). Because of this, our 
more general solution Eq. (27), truncated A^ph = 1 to 
take into account only zero or single microwave photon 
exchange processes, gives identical result to Eq. (39) in 



the rotating frame approach assuming absence of spin- 
flip processes. 

Similarly to Sec. |III A[ we can take into account that 
<C Ep for microwave frequencies which yields pumped 
current in the adiabatic limit (where current is propor- 
tional to uj): 



tRF 



-Tr{rp[a„GT]G«}. 



(40) 



Here all matrices are computed at Ep in the laboratory 
frame after we neglect their frequency dependence in the 
rotating frame by invoking the adiabatic condition uj 
0. 




30 60 90 120 150 180 
Precession Cone Angle (deg) 

FIG. 6: (Color online) The comparison of the dc pumping 
voltage in a clean N|F|I|F|N junction with finite thickness F 
layers (dp = 50) and a clean F|I|F junction whose F layers are 
semi-infinite.^^ The two curves can be computed using ei ther 
the adiabatic NEGF formula in the rotating frame Eq. (40) 
or the adiabatic scattering formula Eq. (43). The parameters 
of these junction are chosen as Ef = —27, A = 27, and 

Ub - 97. 



This expression is equivalent to Eq. (40) due to the fact 



that frequency dependence of all NEGF quantities in the 
latter has been neglected. 



C. Comparison with adiabatic scattering theory 

In the adiabatic limit, one can also employ the Brouwer 
scattering formula^^ which gives the following expression 
for pumped charge current in terms of the derivatives of 
the instantaneous scattering matrix of the device: 



Nevertheless, the application of Eq. (41 ) to clean MTJs 



rST 



Y dtTT l^S,,{EF^t)i^^Sl{E, 

n V <7 



(41) 

We can recast Eq. ( [41] ) in terms of NEGFs for stationary 
transport (which depend on only one energy argument) 
by using the Fisher-Lee formulsl®' for the scattering ma- 
trix 



I6pq + iVfp • G;jEF,t) • Vfq. (42) 



Here G^^ is the submatrix of G 



[E - H(t) - 

which connects edge monolayer of the device attached to 
lead p to the edge monolayer attached to lead q. The 
pumped current is then expressed a d^^ * ^^ * 



rST 



r 

'JdtTi |rpG'-(^f,t)i^[rG"(EF,t)]|.(43) 



with no spin-flip processes, where the scattering matrix 
was obtained directly by matching the wave functions 
across a simple model of Fe|MgO|Fe junction, has pre- 
dicted^^ three orders of magnitude smaller pumping volt- 
age than the rotat ing fr ame formula Eq. ( [4Q| ) applied to 
the same junction) -*^^ ! ^^ Figure [g] explains that the ori- 
gin of this discrepancy is not the particular formalism 
employed, but the assumed MTJ setup which contains 
semi-infinite F leads in Ref.|l6|and finite thickness F lay- 
ers in Ref. I17|18l Although the exact ratio of the voltage 
signals in these two models depends on the chosen values 
of Ef and A at fixed 17^, the voltage signal in F|I|F MTJ 
model always remains below the one in the N|F|I|F|N 
model. 

Both of these models describe unrealistic MTJs — there 
is no any kind of spin accumulation in the F layers in 
the F|I|F model, where semi-infinite F leads simply serve 
to define the spin-dependent scattering states, or spin 
accumulation persists throughout the finite thickness F 
layer in the N|F|I|F|N model. At the same time, the 
prediction of Ref. lISj is far below experimental valued 
(~ 1 nV predicted versus ^ 1 /iV measured at 2 GHz 
FMR frequency), which points out to the need to take 
into account additional ingredients^^ in the MTJ model. 
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FIG. 7: (Color online) The dc pumping voltage in a clean 
F|I|F MTJ with finite thickness F layers (dp = 50) in the 
absence of any SOCs computed using: (i) the exact solution 
Eq. ( |4Q| ) obtained via the rotating frame app roach or, equiva- 
lently, full time-dependent solution Eq. ( 30 ) with one photon 



processes taken into account A^ph = 1; and (ii) truncated (to 
n = ±1) continued fractions solution to double-time- Fourier- 
transformed NEGF equations which gives pumped charge cur- 
rent via Eq. (46). The shaded area marks the interval of 



precession cone angles ^ < 10° beyond which the continued 
fractions solution is not applicable anymore. 



D. Comparison with continued fractions solution to 
double-time-Fourier-transformed NEGF equations 

The double-time-Fourier-transformed NEGF equa- 



tions in Sec. Ill A have been solved before for spin and 
charge pumping problems in an iterative manner using 
continued fractions. Therefore, in this Section we 
discuss the advantage of our solution from Sec. Ill A over 



continued fractions technique using the same F|I|F clean 
MTJs model (with F layers of finite thickness and no 
SOCs) from Fig. [g] as a test bed. 

In the continued fractions method, one starts from the 
equation of motion for Gl^{E) written as 

+G;_i(^)Vtg;(^), (44) 

Here g^(^) = - Hq - 5]^(£;)]-^ is the Fourier trans- 
form of the retarded GF g^(t — in the absence of 
the pumping potential H'(t) = which, therefore, de- 
pends only on the time difference t — t' and can be 
Fourier transformed to a single energy argument. We 
also use notation g;(E) = g^(E + ncj), G^ = 27r(5(0)G; 
and HI = V|g^a iV + VgLi(3-iV'^ . This equation is 
then solvecPSl^^'^ in an iterative manner using Gq = 

[(g5)-i - s5]-i for_n = 0, g;; = G;;_iVtg>„ for 

n > 1, and = G^_^{Vg'^l3n for n < —1. The co- 
efficients c^rl and /3n are generated through continued 
fractions, OLn{l - V'^g^.^iCXn+i^gn) = 1 and /3n(l - 
Vg;_i/3n_iVtg;) = 1. 



The knowledge of GJ^ all ows o ne to express the pumped 
charge current in lead p a d^" * ^^ * 



rCF 



- ^ y dETv{r,{E)G:,{E) 



xr{E + nuj)Gl{E)}[f{E + nto) - f{E)]{45) 

The summation over n in this formula shows how multi- 
photon exchange processes assist current pumping. This 
expression can be used for non-adiabatic external poten- 
tials while in the adiabatic regime uj ^ and at zero 
temperature the difference of Fermi functions is replaced 
by f{E -h nuj) - f{E) ^ ujS{E - Ep), so that only the 
Fermi level states carry the pumped current. 

Although one can in principle solve continued fractions 
for OLn and (5n to arbitrary order n, this is virtually im- 
possible to execute for sizable 3D devices (such as the 
ones in Fig.jl]) due to the need to compute numerous sub- 
matrices of G'^^_i required to obtain GJ^. Instead, most 
of recent applications^^ of the continued fractions so- 
lution to spin pumping in 2DEGs with the Rashba SOC 
have utilized only a few fractions (|n| < 3). In other 



words, the convergence of the sum over n in Eq. (45) can 



be achieved quickly only for small amplitude of the exter- 
nal potential ||V|| which ensures that higher order 
fractions are negligible. The lowest order n = 0, ±1 ver- 
sion of Eq. (45) simplifies to^^ 



+ 00 



^ j dETT{T^{E)g^{E)\^A^{E)\g\E)} 

— oo 

x[f{E + w)-f{E)] 

+ 00 

L I dETT{T,{E)g^{E)VA_r{E)V^g^{E)} 



:[f{E-oj)-f{E)], 



(46) 



where A,(^) = gUE)T{E + nu;)g^{E). 

Figure [tI shows that Eq. (46) is insufficient to ana- 



lyze pumping by magnetization dynamics in MTJs since 
it fails to reproduce the exact solution for dc pumping 
voltage in junctions with no spin flips given by Eq. (39) 
in the rotating frame approach. Because the strength 
of the pumping potential in Eq. (22) is determined by 
Asin6>, pumping voltage computed from Eq. (46) can 
be valid only at small angles cone angles {0 ^10° in 
Fig. [t] this interval would be somewhat larger in F|N 
multilayers). Even at small cone angles, the prediction 
Ip^ (X sin^ ^ stemming from Eq. (^46), which is in 
accord with the "standard model" Eq. (p), becomes in- 
correct in the presence of SOC where /p vs. turns out 
to be quite different (see Sec. IV). We note that one could 
try to use more general Eq. (45), but this would require 



to compute continued fractions and fBn to high order 
n, unlike our non-perturbative solution Eq. (27) which 
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reproduces the exact result in the rotating frame using 
only n = 0, ±1 in the multiphoton GF in Eq. (19). 



IV. THE EFFECT OF INTERFACIAL RASHBA 
SOC ON THE VOLTAGE SIGNAL OF SPIN 
PUMPING IN FIIIN AND FlIlF JUNCTIONS 



Equation (30) applied to clean F|I|N junctions with in- 
terfacial Rashba SOC allows us to understand how the 
dc pumping voltage can appear in such semi-MTJs at 
the adiabatic level. The Rashba SOC is present at the 
F|I interface (i.e., at the last monolayer of the precessing 
F layer that is in contact with the tunnel barrier I) and 
intrinsically participates in the pumping process. This 
is in contrast to other recent theories^^ of spin pump- 
ing in F|N multilayers where SOC is located away from 
the precessing F layer and, therefore, is not essential to 
understand the pumping effect itself. 

Figure IsFa) demonstrates that dc pumping voltage 



OC 



huj in F|I|N junctions emerges as soon as 



the Rashba SOC is "turned on". This could explain 
signal observed experiment all^^^^ in F|I|N junctions, in 
contrast to previous attempts^ ^ which have predicted 
^pump OC (huj)'^. For comparison, Fig.pFb) shows how the 
presence of strong Rashba SOC directly at the F|I inter- 
face also enhances T/^pump in conventional F|I|F MTJs. 

Figure [sf^d) provides additional insight into the charge 
pumping mechanism where we show that the dc pump- 
ing voltage in F|I|N semi-MTJs requires to include ex- 
change of up to ten microwave photons in order to reach 
its asymptotic value. However, since that asymptotic 
value of Fpump is only about 10% higher than the result 
plotted in Figs, [s]^ a) and [sj^c), where only zero or single 
microwave photon exchange processes are taken into ac- 
count, we employ only this lowest order approximation 
in the rest of the paper since [A^ph = 1 in Eq. (30) is 
computationally much less expensive than A^ph = 10\. 

The unique experimentally testable signature of charge 
pumping in F|I|N semi-MTJs that we predict in Fig.jsj^c) 
is angular dependence of Ipump(^) which changes sign 
and it is, therefore, quite different from the usual 
VpumpiO) for F|I|F MTJs shown in Fig. [6] The charge 
pumping in semi-MTJs with weak interfacial Rashba 
SOC can be obtained analytically using the second-order 
perturbation expansion of GF entering Eq. (41) as the 
version of the Brouwer scattering formula: 



Jo 



m{t) X 



dm{t) 



dt 



(47) 



cos^ X - 2 ^^^^ ^ 



sm e/cosc 



Here x is the angle between the axis around which the 
magnetization precesses and the direction of transport 
[in the case of our device in Fig. [ij^b), cos^ X~ ^ sin^ x = 
— 1/2]. The lengthy explicit expression for oc 7^go 
determining the magnitude of the pumped current is pro- 
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FIG. 8: (Color online) The dc pumping voltage in clean F|I|N 
semi-MTJ [panels (a),(c) and (d) and F|I|F MTJ [panel (b)] 
with finite thickness {dp — 50) F layers and non-zero interfa- 
cial Rashba SOC. The Rashba SOC is located within the last 
monolayer of the precessing F layer that is in contact with the 
tunnel barrier I in (a),(c),(d) [as illustrated in Fig. [TJa)], or 
such edge monolayers are present in the left or both F layers 
[as illustrated in Fig. [ijb)] in panel (b). The data in pan- 
els (a),(b),(c) is computed by considering only one microwave 
photon exchange processes, while in panel (d) we show correc- 
tion to this result when up to 9 microwave photons are taken 
into account in Eq. (|30| applied to F|I|N semi-MTJ. 



vided in Appendix [A} In Fig. [sjc), we assume strong in- 
terfacial Rashba SOC (7rso = O.57) so that Ipump(^) vs. 
6 plotted there deviates from this analytical expression 
V^nmi>{0) (X sin^ l9 COS l9/G(l9). Note that for smah TAMR 
[Trso ^ 0-4 according to Fig. |4jc)] G{0) can be con- 
sidered nearly constant, so that V^^rn^iO) 
according to Eq. (47). 



oc sm 



'cost 



The second-order nature of this process can be illus- 
trated using real space Feynman paths where electron im- 
pinging onto the tunnel barrier is reflected with rotation 
of its spin introduced by the Rashba interface. There- 
fore, it has to travel twice through this monolayer to 
reach the right N lead. This picture i s en coded quantita- 
tively in the expression for in Eq. ( A2 ) which contains 
(X 7^so dependence. We note that the same ex 7^so 
angular dependence has also been predicted^^ for linear- 
response STT in F|I|N semi-MTJs with the interfacial 
Rashba SOC, which is in accord with reciprocal nature 
of STT and spin pumping. That is, observation of one of 
these two effects implies, by Onsager reciprocal relations, 
the existence of the other effect.-^ 
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V. DISORDER AND EXTRINSIC SOC 
EFFECTS ON CHARGE PUMPING IN 
MAGNETIC TUNNEL JUNCTIONS 

In this Section, we analyze how disorder and the cor- 
responding extrinsic SOC affects dc pumping voltage in 
conventional F|I|F MTJs. To isolate their effects only, 
we assume that interfacial intrinsic Rashba SOC studied 
in Sec. IV is absent. When extrinsic SOC is negligible, 
we find that dc pumping voltage plotted in Fig. [9|a) is 
increasing in the quasiballistic transport regime (char- 
acterized by the Fano factor F < 1/3 in Fig. [2| and 
then decreases once the diffusive regime (characterized 
by the Fano factor F = 1/3 in Fig. [2]) is reached. At first 
site, this initial increase of Vpump with increasing dis- 
order is counterintuitive, even though conductance also 
decreases with disorder, since pumped current appears 
to be increasing with W. However, it can be explained 
qualitatively as being due to random electron scattering 
in real space which prolongs the average time an electron 
remains in the left F layer where it can interact with pho- 
tons of the microwave pumping field. Similar enhance- 
ment of pure spin current pumping has bee n not iced in 
the diffusive regime in related device setups.^^^^ 

The same disorder used in Fig. [9ja) is related to the 
extrinsic SOC through Eq. which becomes a rele- 
vant effect if Aeso is renormalized by the band structure 
effects to become stronger than its vacuum value by sev- 
eral orders of magnitude. Unlike the interfacial Rashba 



xlO" 



SOC studied in Sec.[IV|which brings novel effects into the 
pumping mechanism, the extrinsic SOC simply reduces 
the dc pumping voltage in F|I|F junctions, as shown in 
Fig. |9jb). Our unified quantum transport treatment of 
spin pumping [Fig.|9|b)] and spin diffusion (Fig.|3]) shows 
that voltage signal of spin pumping in MTJs is brought to 
negligible value when the ratio of the F layer thickness to 
I/gf is dp/Lsf ^ 5. We emphasize that our fully quantum- 
mechanical treatment of the conduction electrons is nec- 
essary to understand such interplay of spin pumping, 
spin accumulation around interfaces^ and spin diffu- 
sion in MTJs since conventional approaclP developed 
for F|N multilayers, where pumping is treated quantum- 
mechanically while subsequent propagation of spins and 
charges is described semiclassically using phenomenolog- 
ical mean free path and spin-diffusion length, is inappli- 
cable to systems containing tunnel barriers where spin 
accumulation is not well-defined.^ 



VI. CONCLUDING REMARKS 

We have derived an exact and efficient for computa- 
tional implementations solution to the equations of mo- 
tion for the double-time-Fourier-transformed NEGFs in 
the presence of time-periodic external potential. Unlike 
continued fractions solutioiP^ for the same equations, 
which is often applied to problems of spirl^ and charge 
pumping^^ by computing only a finite number of con- 
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FIG. 9: The dc pumping voltage in F|I|F MTJs of finite thick- 
ness F layers {dp — 50) with: (a) static disorder of strength 
W within F layers; (b) static disorder of strength = 87 
ensuring diffusive transport regime (see Fig. |2]) and the ex- 
trinsic SOC of strength Aeso determined by such disorder 
via Eq. ([2]). The tunnel barrier I in both panels contains 
binary alloy disorder 5Uh — 0.67. The spin-diffusion length 
corresponding to the values of Aeso is shown in the inset of 
Fig. [I 



tinned fractions while assuming that the amplitude of 
time-periodic external potential is small, our formulas for 
pumped charge Eq. (27) and spin Eq. ([28| currents in the 
leads of a multiterminal devices can be used for arbitrary 
strength of periodic driving potential (thereby covering 
both perturbative^^ and non-perturbative regimes) or 
frequency (thereby covering both adiabatic and non- 
adiabatic pumping regimes). 

This fully quantum- mechanical treatment of pumping 
processes and subsequent propagation of electrons is ap- 
plied to the problem of charge pumping by precessing 
magnetization in the single F layer of F|I|N semi-MTJ 
or F|I|F conventional MTJ in the presence of intrin- 
sic Rashba SOC at the F|I interface. The non-zero in- 
terfacial Rashba SOC, located^^ within the edge mono- 
layer of the precessing F in contact with the tunnel bar- 
rier I, generates non-zero dc pumping voltage in F|I|N 
semi-MTJ at the adiabatic level (i.e., pumping voltage 
is proportional to the microwave frequency uo). This 
could explain observations of voltage signal with such 
properties in the recent experimental^ on microwave- 
driven F|I|N semi-MTJs where previously formulated 
theorieJ^^^ have found only a very small non-adiabatic 
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(oc cj^) voltage signal. We further predict a unique signa- 
ture of this charge pumping phenomenon — the pumped 
charge current changes sign (/ ex sin^ cos for small 
7rso) as the function of the precession cone angle so 
that measuring the corresponding dc pumping voltage 
(^pump sin^^cos^ for small 7rso) would confirm our 
prediction. 

Besides offering quantitative description of charge and 
spin pumping processes, our solution for NEGFs whose 
two energy arguments are connected by the Floquet the- 
orem describing multiphoton emission and absorption 
processes also provides new physical insights: {i) in the 
absence of SOCs, emission or absorption of one photon 
is sufficient to match the exact solution in the rotating 
frame;^^ (n) in the presence of Rashba SOC, exchange 
of up to ten photons is required to reach asymptotic value 
of the pumped currents. Nevertheless, this asymptotic 
value is only about 10% larger than the value obtained 
using just one photon processes in the presence of SOCs. 

We also find that static disorder can increase the dc 
pumping voltage in F|I|F MTJs with finite thickness F 
layers in the quasiballistic transport regime where scat- 
tered electrons spend more time within the precessing F 
layer to interact with microwave photons. The extrinsic 
SOC determined by the impurity potential responsible 
for the diffusive transport regime causes spin relaxation 
which ultimately diminishes the pumping voltage in F|I|F 
MTJs to zero when the spin-diffusion length is about five 
times shorter than the thickness of the F layers. 
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Appendix A: Expression for pumped charge current 
in F|I|N junctions with interfacial Rashba SOC 

Here we provide explicit expression for the coefficient 
which determines the strength of pumped charge cur- 
rent in Eq. (47) for F|I|N semi-MTJ with weak interfacial 
Rashba SOC, as discussed in Sec. [TV| In the absence of 
Rashba SOC, the retarded GF of F|I|N semi-MTJ can be 
written as 



GS = go + gim-(T. (Al) 
Starting from this expression, the second-order perturba- 
tion theory in the powers of 7rso applied to GF in the 
Brouwer scattering formula, recast as Eq. (41), yields 



Z)^ = 4^Re{Tr[rpgoF«giyagirgI 

a 

-T<'g,V^g,Tg\Vagl]}, (A2) 

Here the vector V = (14, F^, 14), which contains the 
strength of the Rashba SOC 7rso, is defined by 
Hrso = • (V X (t) where Hrso is the Rashba Hamil- 
tonian [i.e., the fourth term in Eq. ([T])]. 
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